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Abstract: A key step towards the understanding of the origin of ultra-high energy cosmic 
rays is their mass composition. Primary photons and neutrinos produce markedly different 
showers from nuclei, while showers of different nuclear species are not easy to distinguish. 
To maximise the discrimination with the Pierre Auger Observatory ideally all mass-sensitive 
observables should be combined, but the 10% duty cycle of the fluorescence detector limits 
the use of direct measurements of shower maximum at the highest energies. Therefore, we 
investigate mass-sensitive observables accessible with the surface detectors alone. These are 
the signal risetime in the Cherenkov stations, the curvature of the shower front, the muon-to- 
electromagnetic ratio, and the azimuthal signal asymmetry. Risetime and curvature depend 
mainly on the depth of the shower development in the atmosphere, and thus on primary 
energy and mass. The muon content of a shower depends on the primary energy and the 
number of nucleons, while asymmetry about the shower core is due to geometric effects and 
attenuation, which are dependent on the primary mass. The mass sensitivity of these variables 
is demonstrated and their application for composition studies is discussed. 
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Introduction 

The Pierre Auger Observatory [2, 9] records 
air showers induced by cosmic rays with en- 
ergies > 10 18 eV with unprecedented preci- 
sion and statistics. Key to this are the shape 
of the energy spectrum, the arrival direction 
distribution and the composition of the pri- 
maries. Specifically, a transition of galactic 
to extragalactic cosmic rays could be indicated 
by a marked change of composition (see e.g. 
[8, 5, 4]). While the determination of the di- 
rection and energy of cosmic rays are relatively 
easy , the determination of the identity of pri- 
mary particles from the details of the atmo- 
spheric extensive air showers (EAS) is diffi- 
cult. Differences between neutrino and pho- 
ton induced showers and those initiated by nu- 
clei are relatively easy to detect [3, 1, 6]. The 
differences between nuclear species (e.g. pro- 
ton and iron) are more difficult to measure: on 
average heavier nuclei (of the same total en- 



ergy) develop higher in the atmosphere (lower 
A max ), due to their larger cross-sections and 
the lower energy per nucleon, produce more 
secondary particles, and have a higher muon- 
to-electron ratio at observation level. These 
differences are small and easily concealed by 
the variations due to the primary energy or by 
uncertainties in hadronic interaction models, 
and are covered by shower fluctuations and the 
limited precision of the measurement with a re- 
alistic detector. The Pierre Auger Observatory 
is a hybrid detector, combining a large array of 
water-Cherenkov detectors (surface detector, 
SD, with 100% duty cycle) and 24 fluorescence 
telescopes (FD, with 10% duty cycle). Hybrid 
events have a good energy estimate and pro- 
vide a direct measurement of X max . But 90% 
of the events have information only from the 
surface array and, therefore, observables from 
the array are needed for composition analysis 
at the highest energies. In the following we 
discuss four observables which all rely on the 
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time structure of the shower front as measured 
with a 25 ns flash ADC system (FADC) of the 
water-Cherenkov detectors. They are all sen- 
sitive to the primary mass, via their depen- 
dence on the height of shower development and 
the muon content in the shower: (i) The sig- 
nal risetime in water tanks at some distance 
from the shower core depends mainly on the 
height of the shower development: the higher 
the shower development in the atmosphere, the 
smaller the time spread between the particles 
and the smaller the risetime. (ii) The cur- 
vature of the shower front reflects the height 
of the shower development: the further the 
main shower development from the detector, 
the flatter the shower front, (iii) The muon 
content of a shower can be assessed by looking 
for tell tale signs of sudden, large energy de- 
posits in the time traces, as they are expected 
from muons: the frequency of large jumps in 
the signal from one FADC bin to the next re- 
lates to the abundance of muons. (iv) The 
azimuthal asymmetry in the signal risetimes 
within the shower plane is an effect of a vary- 
ing muon-to-electromagnetic component, and 
as such it depends on the shower development 
and the muon content. 

While all of these variables could in principle 
be used for composition studies with a perfect 
detector, there are limitations due to the finite 
accuracy with which they can be measured. 
Therefore, it needs to be shown that measure- 
ment uncertainties are small enough that the 
mass differences can be detected reliably. Once 
the mass sensitivity of the individual variables 
is established, a multivariate analysis method 
will be needed to combine all mass-sensitive 
observables to maximise the mass separation 
power. A combination of risetime and curva- 
ture has already been used successfully to pose 
an upper limit on primary photons [1] . 

The risetime of the signal, t 1 / 2 

The signals in the 10 m 2 water-Cherenkov de- 
tectors are characterised by fi/2, the time to 
rise from 10 to 50% of the integrated signal. 
The signals, up to 3 fis long, are recorded with 



a 25 ns FADC system. The risetime in an indi- 
vidual station depends on core distance, zenith 
angle and energy and has been parameterised 
from experimental data to correct events from 
different zenith angles. The measurement un- 
certainty otj/j as function of signal S in the 
tank, the core distance r and zenith angle 9 
has been determined from a small number of 
twin tanks and from pairs of detectors at sim- 
ilar core distances (within 100 m). A param- 
eterisation of c tl/2 is used to define the error 
bars shown in fig. 1 (top). Measurements of 
tx/2 for signals < 20 VEM are severely influ- 
enced by Poisson fluctuations due to low num- 
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Figure 1: Top: (A) derived from a single 
event. The black line is the predicted rise- 
time while the data points represent measure- 
ments with the error bars from twin station 
studies. Combining each stations' A from the 
expectation an average for the event can be de- 
fined. Bottom: (A) as a function of A max , for 
nearly-vertical hybrid events. 
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bers of particles and are not used. A correction 
of the risetime for the azimuthal asymmetry in 
the shower, particularly important between 35 
and 50°, is applied. For each station the de- 
viation A of the measured risetime, in units 
of the accuracy of measurement, from the av- 
erage values found for r and 9 at a fixed ref- 
erence energy E « 10 19 eV is formed, thus: 

<rt 1/2 (S,r,0) 

event are then combined in an average event 
deviation (A) as seen in fig. 1 (top), which 
should be larger for showers developing deeper 
than the reference and smaller for those devel- 
oping higher. Using hybrid events [10] it can 
be tested whether the so-constructed variable 
(A) indeed correlates with X max . The data 
have been binned in X max , and for all events in 
one A max bin the average of (A) (i.e. (A)) has 
been formed. This is shown in fig. 1 (bottom). 
A linear dependence is observed which allows 
estimation of A max from events observed with 
the SD alone. 
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Shower Front Curvature 

The trigger times of the array detectors away 
from the core lag behind what one would ex- 
pect from a plane of particles traveling at the 
speed of light, indicating the shower front is 
curved. With the assumption of a hemispheri- 
cal shower front a radius of curvature, R c , can 
be derived. The radius of curvature is an SD 
observable that Monte Carlo simulations show 
is mass sensitive (see Fig. 2). Larger dis- 
tances of A max from the array yield larger radii 
of curvature, while deeply penetrating show- 
ers produce shower fronts with smaller radii 
of curvature. Also muons tend to travel in 
straight lines and electrons are more scattered, 
so muons arrive at the detectors first, and 
with a narrower time spread, as compared to 
the electro-magnetic component. Therefore, a 
shower with a larger number of muons will pro- 
duce station triggers that are compact in time, 
resulting in a flatter shower front than for a 
similar shower that is muon-poor. The radius 
of curvature is strongly dependent on zenith 
angle because the distance from the array to 



Figure 2: Top: Comparative radii of curvature 
from simulations of three primaries (solid-blue: 
proton, solid-red: iron, dashed-magenta: pho- 
ton). Bottom: R c is shown as a function of 
X — A max for nearly vertical hybrid showers 
(1.1 < sec(0) < 1.2). The line is a linear fit 
to the data points correlating R c with the dis- 
tance from the detector to the shower maxi- 
mum. 

shower maximum increases with zenith angle 
(oc (Xo sec (9 — A max )). Fig. 2 demonstrates the 
typical resolution between the radius of curva- 
ture and changes in X — A max . 

Muon Content 

As a typical muon deposits much more energy 
(w 240 MeV) in a water tank than an electron 
or photon (w 10 MeV), spikes are produced 
over the smoother electromagnetic background 
in the FADC time traces. Thus, muons mani- 
fest themselves as sudden variations in the sig- 
nal AV from one FADC bin to the next. The 
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Figure 3: Distribution of jumps for Auger 
events with energy between 8-12 EeV and 40- 
50°, averaged over 251 stations - 100 showers - 
within the same core distance interval. 

expected distributions of AV for purely elec- 
tromagnetic and muonic traces look very dif- 
ferent. A mix of electromagnetic and muonic 
signals is needed to fit the measured distribu- 
tion (see Fig. 3). The number of muons in 
the shower is then estimated via the number 
of AVs above 0.5 VEM with = aN(AV > 
0.5VEM); with a w 1.4 being only weakly de- 
pendent on core distance and zenith angle. The 
resolution is about 25%, systematic errors 
are below 10%. 

Azimuthal Asymmetry 

The observed azimuthal asymmetry in the rise- 
times offers another handle on primary compo- 
sition determination. Its magnitude is a mea- 
sure of the degree of shower development and 
hence it is sensitive to primary composition 
[7]. The dependence of the risetimes with az- 
imuthal angle £ is fitted using ti/ 2 (r, £) = a + 
b cos £, for all stations with 500 m < r < 2000 
m and with signals greater than 10 VEM. The 
asymmetry factor - is sensitive to the distance 
between the detector and the shower maximum 
(shown in fig. 4 as sec 9 oc X — X max ). Max- 
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Figure 4: Development of the asymmetry fac- 
tor in the risetime with zenith angle. The pri- 
mary energy is ~ 10 18 eV. 

imal asymmetry occurs at a unique value of 
X — X max found most frequently at a corre- 
sponding zenith angle. Monte Carlo simula- 
tions show that this zenith angle (0 asyrnax ) is 
different for proton and iron showers (Fig. 4) 
reflecting the systematic difference in the aver- 
age A max . 

References 

[1] J. Abraham et al. [Pierre Auger Collabo- 
ration]. Submitted to Astrop. Phys. 

[2] J. Abraham et al. [Pierre Auger Collabo- 
ration]. Nucl. Instr. Meth. A, 523:50-95, 
2004. 

[3] J. Abraham et al. [Pierre Auger Collabo- 
ration]. Astrop. Phys., 27:155, 2007. 

[4] D. Allard et al. Astrop. Phys., 27:61, 2007. 

[5] V. Berezinsky et al. Phys. Rev. £>., 74:61, 
2006. 

[6] O. Blanch [Pierre Auger Collaboration]. 

In these proceedings (#603). 
[7] M. T. Dova [Pierre Auger Collaboration]. 

In Proc. 29 th ICRC (Pune), pages 101- 

104, 2005. 

[8] A. M. HiUas. J. Phys. G, 31:95, 2005. 

[9] T. Suomijarvi [Pierre Auger Collabora- 
tion]. In these proceedings (#299). 
[10] M. Unger [Pierre Auger Collaboration]. In 
these proceedings (#594). 



